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Abstract 

Magnetic relaxation and critical current density have been measured on a 
MgB2 thin him in a wide region of temperature with the magnetic field up 
to 8 T. The irreversibility line has also been determined. It is found that the 
relaxation rate has a very weak temperature dependence below 1 / 2Tc showing 
a clear residual relaxation rate at zero K, which cannot be easily explained as 
due to thermally activated fiux creep. Furthermore the relaxation rate has a 
strong field dependence. The fiux dynamics of thin films are very similar to 
that of high pressure synthesized bulks although the relaxation rate in thin 
film is systematically higher than that of a bulk sample. All the results here 
together with those from the bulk samples suggest that the flux dynamics 
may be dominated by the quantum effects, such as quantum fluctuation and 
tunneling. 

74.25.Bt, 74.20.Mn, 74.40.+k, 74.60.Ge 



Typeset using REVTgX 



I. INTRODUCTION 



The recently discovered new superconductor MgB2 becomes very attractive due to its 
potential applications HTj-Q- One important issue concerns however in which region on the 
field-temperature ( H-T ) phase diagram it can carry large critical current density ( jc ) 
and thus can be used in the future for industry. This jc is controlled by the mobility of the 
magnetic vortices, and vanishes above the melting line between the vortex solid and liquid. 
This melting can be induced by strong fiuctuation of the vortex position by either thermal 
effect or quantum effect. At T = K only the quantum fiuctuation is left. A finite linear 
resistivity pun = {E/j)j^o will appear above this melting point showing the reversible fiux 
motion. Most of the published results on the fiux dynamics of MgB2 system were obtained 
from bulk samples l^j-Q . In this paper we present the investigation on the fiux dynamics on 
an MgB2 thin film by dynamical magnetic relaxation method. 



II. EXPERIMENTAL 



The thin films of MgB2 were fabricated on (1102) AI2O3 substrates by using pulsed laser 
deposition technique |TD[. The amorphous B thin film was first deposited and then it was 



sintered at a high temperature in Mg vapor. The films are typically 400 nm thick with 
predominant c-axis orientation ( c-axis is perpendicular to the film surface ). A rectangular 
sample of size 2.1mm x 4.9mm was chosen for the magnetic sweeping measurement. The 
magnetic measurements were carried out by a Quantum Design superconducting quantum 
interference device ( SQUID, MPMS 5.5 T ) and a vibrating sample magnetometer ( VSM 
8T, Oxford 3001 ) at temperatures ranging from 2 K to Tc and external field up to 8 T. 
For the magnetic sweeping measurement the M(H) curve was measured with different field 
sweeping rates ( 0.005 T/s to 0.02 T/s ) and integral time of 60 ms to 960 ms. The pressure 
of He gas in the sample chamber for thermal exchange was kept at 0.4 bar during the 
measurement. The superconducting transition is sharp with the transition temperature Tc 



2 



of about 38 K as observed from the temperature dependence of magnetization. 



III. RESULTS 

In Fig.l we show the magnetization hysteresis loops ( MHL ) measured at temperatures 
ranging from 2 K to 37 K. The symmetric MHLs observed at temperatures up to 37 K 
indicate the dominance of the bulk current instead of the surface shielding current. The 
MHLs measured at low temperatures, such as 2 K to about 10 K show dense and small flux 
jumps in the low field region leading to a smearing of the superconducting critical current 
density. It is this effect that makes the width of the magnetization AM below 10 K be 



anomalously smaller than that at 14 K, which is discussed elsewhere [|rT|. From these MHLs 
one can calculate jc via jc = 20AM/V"a(l — a/36) based on the Bean critical state model, 
where V, a and b are the volume, width and length of the sample, respectively. The result 
of jc is shown in Fig. 2. It is clear that the magnetic critical current density jc of our sample 
is rather high. For example, at T = 18 K and fiQH = 1 T, we have jc = 2 x lO^A/cm^, 
which is about one order of magnitude higher than that of high pressure synthesized bulk 
samples i,0,|ll. 

For investigating the flux dynamics the jc{H) curves have been measured with two dif- 
ferent field sweeping rate 0.02 T / s and 0.01 T / s. It is interesting to note that although 
for bulk samples [H,^ there is a small tail on each MHL in high field region, here on the 
thin film sample there is no such tail. Accordingly the jc{H) curves show neither a small 
tail in high field region. This further corroborates our earlier suggestion that the small tail 
observed in bulk samples is due to some secondary effect, such as some local regions with 
very strong pinning or the surface pinning by the tiny grains. In the present thin film sample 
with much better uniformity this effect is certainly disappeared. From the contour of jc vs. 
H shown in Fig. 2 one can see that all the curves have two different regions. In the low field 
region jc drops slowly with H. When the field is increased further and exceeded a threshold 
the jc decreases drastically showing the gradual setting in of the reversible flux motion. One 
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can determine the phase transition hne which separates the irreversible and the reversible 
flux motion by taking a criterion, here for example jc = 1000 A/cm'^. The same criterion 
was used by Bugoslavsky et al. for bulk samples 0]. Since our measurement was done at a 
maximum field of 8 T, for low temperatures we use a reasonable extrapolating way, i.e., to 
follow the tendency of jc vs. H at a higher temperature, for example at 10 K or 14 K, down 
to the criterion jc = 1000 A/cm'^ to derive the irreversibility field. This method, although 
with some uncertainties for the values of Hirr{T) at low temperatures, will lead us to derive 
an almost complete curve of Hirr{T). The error bar of Hirr{T) at 2 K is about ± 0.5 T. 
The irreversibility lines by following this method is shown in Fig. 3 together with that from a 
high pressure synthesized bulk sample. Similar to what found in bulk samples , it 

is easy to see in the present thin film that the Hirr{T) extrapolates to a rather low field, for 
example, Hirr{0) ~ 9.2± 0.5 T, while the Hc2{T) extrapolates to a much higher value ( ~ 
15 T ) |13[ at zero K. There is a large separation between the two fields Hc2{0) and Hirr{0). 
This effect observed both in bulk and thin film samples may suggest that the relatively low 
Hirr{T) in MgB2 is not due to the easy flux motion through some weak pinning channels, 
rather it reflects probably a more intrinsic property, specialy in a rather clean system. 



IV. DISCUSSION 

A. Large separation between Hirr{0) and Hc2{0) and possible evidence for the 

quantum vortex liquid in MgB2 

If following the hypothesis of the vortex liquid above Hirr{T), we would conclude that 
there is a large region for the existence of a quantum vortex liquid at zero K. This can be 
attributed to a quantum fluctuation effect of vortices in bulk MgB2- Although the lowest 
temperature in our present experiment is 2 K, however, from the experimental data one 
cannot find any tendency for Hi^j.{T) to turn upward to meet the -ffc2(0) at zero K. This may 
indicate that the vortex melting in our present film is due to the strong quantum fluctuation 



4 



which smears the perfect vortex lattice leading to the vanishing of the shear module Cqq 
of the vortex matter ( probably within grains ). Dense disorders will strengthen the shear 
module and thus enhance the irreversibility line, that's why the ifj^r in present thin film is 
higher than that in bulks ( shown in Fig.3 ). The irradiation of protons by Bugoslavsky et 



al. [Ill] did not suppress but strongly increase jc at a high field would suggest that the low 
value of Hirr (T) measured in unirradiated bulks and present thin film is not due to the weak 
links since otherwise the jc value would drop even faster with increasing the magnetic field 
after the irradiation. 

In order to investigate the flux dynamics in the vortex solid state below Hirr we have 
carried out the dynamical magnetic relaxation measurement. Assuming an uniform current 
density over the cross-section of a superconducting ring, one can determine the supercon- 
ducting current density jc from the magnetic moment M via 

M = ^7TJcd{Ro^ - R,') (1) 

where d is the film thickness, Ro and Ri are the outer and inner radius of the ring, respec- 
tively. It is important to note that the magnetic moment M in Eq.(l) is understood as being 
due to only the superconducting current excluding any additional contribution from an equi- 
librium magnetization. In other words, M is obtained by subtracting from the ZFC magnetic 
moment with the equilibrium magnetic moment as determined from the FC process, i.e., 
M = M ( measured ) - M ( equilibrium ). For a superconducting thin film the equilibrium 
magnetic moment is normally negligible since the volume of the superconducting material 
is very small. Based on different external conditions, there are two techniques to measure 
the magnetic relaxation, namely the so-called conventional and dynamical relaxation. The 
so-called conventional relaxation is to measure the time dependence of the superconducting 
current density jc at a certain temperature and field [p!5|-p!7[]. After a waiting time when 
the magnetic field is fixed ( or say field sweeping is stopped ), the first data point is taken. 
For a conventional relaxation measurement, total observation time should be very long. The 
second method is the so-called dynamical relaxation, i.e., to measure the MHL with different 



field sweeping rate. One can easily understand the difference between these two methods 
from the electromagnetic response of a superconducting ring. The electromotive force in a 
ring is 

Ex2nR = nR' x - wdL^ (2) 

dt dt ^ ^ 

where E is the electric field established within the ring, R = {Ro+Ri)/2, H is the external 
field, L = fioR[ln{8R/w) — 1/2] is the self-inductance of the ring flSl, w{= Ro — Ri) is the 



width of the ring. For the relaxation process, since dH/dt = 0, the electric field can be 
determined by 

E = -f^[lni^-^)-^^ (3) 
2n ^ w ' 2^ dt ^ ^ 

This method is inapplicable when the irreversible magnetic signal is comparable to the 
equilibrium magnetization. For MgB2 with a very narrow MHL this method shows a clear 



drawback. However one can choose the so-called dynamical relaxation method ||T9| , pO| . This 
technique can be accomplished by using the sensitive VSM or a torque magnetometer. In a 
field sweeping process, if the field sweeping rate is high enough, the last term in eq.(2) can 
be neglected, therefore the electric field E can be determined by 

RdijioH) 

and jc can be determined from the width of the micro-hysteresis loops around a certain 
field via jc = 20AM/Va{l — a/36) based on the assumption that the current density jc 



is uniform throughout the cross section of the ring Since the vortices are forced to 

move by the external field sweeping, this process is thus called as the dynamical relaxation 
||T9| , |20[] . For a superconducting disc the magnetic moment is contributed mainly from the 
current circulating near the perimeter of the sample. For example, the 7/8 of the total 
magnetic moment is contributed by the current flowing from 1/2R to R of a superconducting 
disc, where R is the radius of the disc. Therefore from a rough estimation, it is safe to 



derive the E(j) or V(I) relation for a disc sample by using Eq.(4). As indicated in ||20 



and [0, the normalized relaxation rate can be determined via Q = dlnjs/dlnE in the 



dynamical relaxation process, and it should be identical to 5 = —dlnM / dint determined 



in the conventional relaxation process. Later on Wen et al. p2|, Perkins and Caplin |23 



and Ji et al. ||2^ have shown that the conventional relaxation, dynamical relaxation and the 
DC transport method should give the same information of flux motion although the voltage 
range is different. ( voltage in the dynamical relaxation is much higher than that in the 
conventional relaxation ). 

The raw data with two different field sweeping rate ( 0.02 T / s, 0.01 T / s) are shown 
in Fig. 2. The Q values vs. field for different temperatures are determined and shown in 
Fig. 4. It is clear that the relaxation rate increases monotonically with the external magnetic 
field and extrapolates to 100% at about the melting point H.i„- At 2 K it is found from 
the Q(H) data that the melting field ( where Q = 1 ) is about 8.7 ± 0.5 T, being close to 
Hirr{T = 2K) ^ 9.2T ± 0.5 determined from the jc{H) curve. It is known that the Hirr{T) 
is rather stable in low temperature region, therefore we can anticipate a rather low value of 
Hirr{0) which is below 10 T being much lower than Hc2{0). As already pointed out in our 
earlier publication the large separation between H^^J''{0) and -ffc2(0) may manifest the 
existence of the quantum vortex liquid due to strong quantum fluctuation of vortices in the 
pure system of MgB2- 

Theoretically, quantum melting of the vortex solid has been proposed by some authors 
and preliminarily verified by experiments. [^,^. Solid evidence is, however, still 



lacking mainly because either the values of ifjrr(O) and Hq{T) are too high to be accessible, 
such as in the classical Chevrel phase PbMoS system [^, or the separation between them 



is too small [28,29| leading to a difficulty in drawing any unambiguous conclusions. Here we 



try to have an rough consideration on the quantum melting field proposed by Blatter 



et al. |^5|] for 2D system, 

HU0)/H,2{0) = 1 - l.2exp{-n^ClRQ/AR2D) (5) 
where Cl is the Lindermann number, Rq = ^/e^ ~ A.lkQ, R2D is the sheet resistance. 



Since the new MgB2 sample has a much higher charge density and thus a much lower sheet 
resistivity, according to above relation, Hm should be more close to i/c2(0) comparing to 
high — Tc superconductors ( HTS ). This is in contrast to the experimental observations 
which may be explained as that the MgB2 is not a quasi-2D system. Another approach was 
proposed by Rozhkov and Stroud PT| , 

H^{0)/H,2{0) = Bo/ {Bo + ife2(0)) (6) 

with Bo = /?mpC^s$o/47rA(0)^g^, where s is the spacing between layers, rrip is the pair 
mass, q the pair charge ( ~ 2e ), C the light velocity, A(0) the penetration depth at zero 
K, /5 ~ 0.1. If comparing again the present new superconductor MgB2 with HTS, A(0), q 
and f3 are more or less in the same scale, the difference comes from rup and s. Therefore a 
preliminary conclusion would be that in MgB2 either the pair mass rUp or the layer spacing 
s is much smaller than that of HTS. 

B. Residual relaxation rate at zero K and weak temperature dependence of the 

relaxation rate 

A strong quantum fluctuation normally favors a strong quantum tunneling creep. In 
order to see that, we plot in Fig. 5 the temperature dependence of the relaxation rate Q. It 
is clear that there is a clear residual relaxation rate for all fields and the relaxation rate in 
wide temperature region keeps rather stable against the thermal activation and fluctuation. 
According to the thermally activated flux motion model 

E = voBexpi-^^^^) (7) 

where E is the electric field due to TAFM over the activation energy U{jc,T), vq is the 
average velocity of the flux motion, B is the magnetic induction. The relaxation rate is 

Q = ^ = _i.,r(M%Ii^) (8) 
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For any kind of U{jc) relation, a finite slope of dU/djc is expected. Therefore a much 
stronger temperature dependence of Q should be expected for thermally activated flux mo- 
tion. This is in contrast to the experimental data. However, when the melting point Hirr{T) 
is approached the relaxation rate will quickly jump to 100%. This may indicate that the 
thermal fluctuation is not the dominant process for the flux depinning in the new super- 
conductor MgB2. It shows high possibility for the vortex quantum melting even at a finite 
temperature. Worthy to note is that the quantum tunneling rate is extremely low at an 
intermediate field, such as Q = 0.3% at 2 K and 2 T, but is rather high at a high field, for 
example Q = 20% at 2 K and 7 T. This may imply that the field will greatly enhance the 
vortex quantum fluctuation and tunneling. 

The small relaxation rate at a relatively low field has also been measured by Thompson 
et al. [|] who regarded it as a highly stable superconducting current density in MgB2. 
Actually the relaxation rate can be rather high when the magnetic field is increased to a 
higher value. The extremely small relaxation rate and weak temperature dependence at a 
finite temperature at a low field is probably induced by a strong pinning barrier relative 
to the thermal energy, i.e., ksT << Uc, where Uc is the intrinsic pinning energy. Recently 
it was concluded P^ , |55[] that the Uc is in the scale of 1000 K being much higher than 
the thermal energy fc^T. Therefore for the new superconductor MgB2 the pinning well 
is too deep leading to a trivial influence of the thermal activation and fluctuation. It thus 
naturally suggests that the quantum fluctuation and tunneling plays an more important role. 
Therefore, together with the fact discussed in last subsection, it is tempting to suggest that 
at a finite temperature the melting between a vortex solid and a liquid is due to quantum 
fluctuation instead of the thermal fluctuation. 



V. CONCLUSION 

In conclusion, in a thin film sample of MgB2, the flux dynamics and the irreversibility 
field are investigated. Just like in the bulk samples, is is found that the irreversibihty 
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field is rather low comparing to the upper critical field in low temperature region showing 
the possible existence of the quantum vortex liquid due to strong quantum fiuctuation. 
The weak temperature dependence but strong field dependence of the relaxation rate may 
further suggest that the vortex melting at a finite temperature is also induced by the strong 
quantum fiuctuation. The reason for such a strong quantum effect is still unknown, but 
it may be related to the superconducting mechanism of MgB2, such as the relatively low 
upper critical field. 



ACKNOWLEDGMENTS 

This work is supported by the National Science Foundation of China (NSFC 19825111) 
and the Ministry of Science and Technology of China ( project: |nkbrsf-gl999064| 602 ). HHW 
gratefully acknowledges Prof. B. Ivlev and Dr. A. F. Th. Hoekstra for fruitful discussions, 
and continuing financial support from the Alexander von Humboldt foundation, Germany. 



10 



[1 
[2: 

[4] 
[5] 

[6' 

[7 



[9 

[10 
[11 



REFERENCES 

E-mail address: hhwen@aphy.iphy.ac.cn. 

J. Nagamatsu, N. Nakagawa, T. Maranaka, Y. Zenitani and J. Akimitsu, Nature 410, 
63(2001). 

Y. Bugoslavsky, G. K. Perkins, X. Qi, L. F. Cohen, and A. D. Caphn, Nature 410, 
563(2001). 

D. K. Finnemore, J. E. Ostenson, S. L. Bud'ko, G. Lapertot, and P. C. Canfield, Cond- 
mat / 0102114. 

D. C. Larbalestier, et al.. Nature 410, 186(2001). 

J. R. Thompson, M. Paranthaman, D. K. Christen, K. D. Sorge, H. J. Kim and J. G. 

Ossandon, Cond-mat / 0103514 ( 2001 ). 

H. H. Wen, S. L. Li, Z. W. Zhao, Z. A. Ren, G. C. Che and Z. X. Zhao, Chin. Phys. 
Lett. 86, 816 ( 2001 ). 

H. H. Wen, S. L. Li, Z. W. Zhao, Z. A. Ren, G. C. Che and Z. X. Zhao, submitted to 
Physica C. 

M. J. Qin, X. L. Wang, H. K. Liu, S. X. Dou, Cond-mat / 0104112 ( 2001 ). 

M. Dhalle, P. Toulemonde, C. Beneduce, N. Musohno, M. Decroux, Cond-mat / 0104395 
( 2001 ). 



W. N. Kang et a/.. Science in press, |cond-mat/0103179| (2001) 



Z. W. Zhao, H. H. Wen, S. L. Li, Y. M. Ni,H. P. Yang, W. N. Kang, H. J. Kim, E. M. 
Choi and S. L Lee, Cond-mat / 0104249 ( 2001 ). 



[12] Y. Takano, H. Takeya, H. Fujii, H. Kumakura, T. Hatano, K. Togano, H. Kito, and H. 
Ihara, Cond-mat / 0102167. 



11 



[13] S. L. Bud'ko, C. Petrovic, G. Lapertot, C. E. Cunningham, and P. C. Canfield, Cond- 
mat / 0102413. 



[14 

[is: 

[16 



[ir; 

[18 
[19 
[20 

[21 

[22 

[23 

[24 
[25 

[26 
[27; 



Y. Bugoslavsky et al, |cond-mat/010415q (2001). 



Y. Yeshurun, A. P. Malozemoff, A. Shaulov, Rev. Mod. Phys. 68, 911 (1996). 

J. R. Thompson, Yang Ren Sun, L. Civale, A. P. Malozemoff, M. W. McElfresh, A. D. 
Marwick and F. Holtzberg, Phys.Rev.B 47, 14440(1993). J. R. Thompson, Yang Ren 
Sun, D. K. Christen, Phys. Rev. B 49, 13287(1994). 

P. J. Kung, M. P. Maley, M. E. McHenry, and J. O. WiUis, M. Murakami, and S.Tanaka, 
Phys. Rev. B 48, 13922(1993). 

J. Gilchrist and E. H. Brandt, Phys. Rev. B 54, 3530(1996). 

M. Jirsa, L. Pust, H. G. Schnack, R. Griessen, Physica C 207, 85(1993). 

H. G. Schnack, R. Griessen, J. G. Lensink, C. J. van der Beek, and P. H. Kes, Physica 
C 197, 337(1992). 

C. J. van der Beek, G. J. Nieuwenhuys, P. H. Kes, H. G. Schnack, and R. Griessen, 
Physica C 197, 320(1992). 

H. H. Wen, H. G. Schnack, R. Griessen, B. Dam, and J. Rector, Physica C 241, 
353(1995). 

G. K. Perkins, A. D. Caplin, Phys. Rev. B 54, 1255(1996). 

H. L. Ji, Z. X. Shi, Z. Y. Zeng, X. Jin, X. X. Yao, Physica C217, 127(1993). 

G. Blatter, and B. Ivlev, Phys. Rev. Lett. 70, 2621 ( 1993 ). G. Blatter, and B. I. Ivlev, 
Phys. Rev. B 50, 10272 ( 1994 ). G. Blatter, et al. Phys. Rev. B 50, 13013 ( 1994 ). 

A. Kramer, and S. Doniach, Phys. Rev. Lett. 81, 3523 ( 1998 ). 

R. Ikeda, Int. J. Mod. Phys. B 10, 601 ( 1996 ). 



12 



[28] T. Sasaki, W. Biberacher, K. Neumaier, W. Hehn, A. Andres, Phys. Rev. B 57, 10889 
( 1998 ). 

[29] S. Okuma, Y. Imamoto, and M. Morita, Phys. Rev. Lett.86, 3136(2001). 

[30] C. Rossel, E. Sandvold, M. Sergent, R. Chevrel, and M. Potel, Physica C 165, 233(1990). 
C. Rossel, O. Pena, H. Schmitt, and M. Sergent, Physica C181, 363 ( 1991 ). 

[31] A. Rozhkov, and D. Stroud, Phys. Rev. B 54, R12697 ( 1996 ). 

[32] Z. W. Zhao, H. H. Wen, S. L. Li, Y. M. Ni, Z. A. Ren, G. C. Che, H. P. Yang, Z. Y. 
Liu and Z. X. Zhao, Chinese Physics 10, 340(2001). 

[33] H. Jin, H. H. Wen, S. L. Li, Z. W. Zhao, Y. M. Ni, Z. A. Ren, G. C. Che, H. P. Yang, 
Z. Y. Liu and Z. X. Zhao, Chin. Phys. Lett. 18, 823 ( 2001 ). 



13 



FIGURES 

FIG. 1. Magnetization hysteresis loops measured at 2, 4, 6, 8, 10, 14, 18, 22, 26, 30, 32, 35, 
37, K and 38 K ( from outer to inner ). All curves here show a symmetric behavior indicating 
the importance of bulk current instead of surface shielding current. The MHLs measured at low 
temperatures ( e.g., 2 K to 10 K ) are too close to be distinguishable. Dense and small flux jumps 
have been observed below 10 K near the central peak. 

FIG. 2. Critical current density jc calculated based on the Bean critical state model. At each 
temperature the data has been measured with two field sweeping rate: 0.02 T / s and 0.01 T / 
s. The faster sweeping rate corresponds to a higher dissipation and thus higher current density. 
From these data one can calculate the dynamical magnetic relaxation rate Q. The jc{H) curves 
measured at low temperatures are very close to each other showing a rather stable value of Hi^r 
when T approaches zero K. We use a criterion of jc = 1000^/cm^ to determine the irreversibility 
line. 

FIG. 3. H-T phase diagram for the new superconductor MgB2. The circles represent the 

bulk irreversibility lines Hi„. of bulk samples: open and filled circles represent two measurement 
by VSM on two bulk samples. The open diamond symbols represent the irreversibility line of the 
thin film. The squares represent the upper critical field Hc2{T): filled squares are from resistive 
measurement; open squares are from the M(T) measurement by SQUID. All lines are guide to the 
eye. 

FIG. 4. Field dependence of the relaxation rate at temperatures of 2, 4, 6, 8, 10, 14, 18, 20, 
22, 24, 26, 28, 30, 32, 35 K. The dashed line is a guide to the eye for 2 K. It is clear that Q wih 
rise to 100% at about 8.7 T at 2 K. Since Hirr is rather stable at low temperatures, it is safe to 
anticipate that Hirr{0) < lOT being much smaller than -ffc2(0) « 15 T. 
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FIG. 5. Temperature dependence of the relaxation rate at fields of 1 to 8 T with increments 
of 1 T. A clear residual relaxation rate is observed at all magnetic fields. The relaxation is weakly 
dependent on temperature until the irreversibility temperature is reached. These are difficult to 
be understood in a frame of thermal acitivation flux motion and thermal depping. 
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